A macroevolutionary role for chromosomal fusion and fission in Erebia butterflies.

The impact of large-scale chromosomal rearrangements, such as fusions and fissions, on speciation is a long-standing conundrum. We assessed whether bursts of change in chromosome numbers resulting from chromosomal fusion or fission are related to increased speciation rates in Erebia, one of the most species-rich and karyotypically variable butterfly groups. We established a genome-based phylogeny and used state-dependent birth-death models to infer trajectories of karyotype evolution. We demonstrated that rates of anagenetic chromosomal changes (i.e., along phylogenetic branches) exceed cladogenetic changes (i.e., at speciation events), but, when cladogenetic changes occur, they are mostly associated with chromosomal fissions rather than fusions. We found that the relative importance of fusion and fission differs among Erebia clades of different ages and that especially in younger, more karyotypically diverse clades, speciation is more frequently associated with cladogenetic chromosomal changes. Overall, our results imply that chromosomal fusions and fissions have contrasting macroevolutionary roles and that large-scale chromosomal rearrangements are associated with bursts of species diversification.

Repeated upslope biome shifts in Saxifraga during late-Cenozoic climate cooling.

Mountains are among the most biodiverse places on Earth, and plant lineages that inhabit them have some of the highest speciation rates ever recorded. Plant diversity within the alpine zone - the elevation above which trees cannot grow-contributes significantly to overall diversity within mountain systems, but the origins of alpine plant diversity are poorly understood. Here, we quantify the processes that generate alpine plant diversity and their changing dynamics through time in Saxifraga (Saxifragaceae), an angiosperm genus that occurs predominantly in mountain systems. We present a time-calibrated molecular phylogenetic tree for the genus that is inferred from 329 low-copy nuclear loci and incorporates 73% (407) of known species. We show that upslope biome shifts into the alpine zone are considerably more prevalent than dispersal of alpine specialists between regions, and that the rate of upslope biome shifts increased markedly in the last 5 Myr, a timeframe concordant with a cooling and fluctuating climate that is likely to have increased the extent of the alpine zone. Furthermore, alpine zone specialists have lower speciation rates than generalists that occur inside and outside the alpine zone, and major speciation rate increases within Saxifraga significantly pre-date increased rates of upslope biome shifts. Specialisation to the alpine zone is not therefore associated with speciation rate increases. Taken together, this study presents a quantified and broad scale perspective of processes underpinning alpine plant diversity.

Phylogenomic analyses of Sapindales support new family relationships, rapid Mid-Cretaceous Hothouse diversification, and heterogeneous histories of gene duplication.

Sapindales is an angiosperm order of high economic and ecological value comprising nine families, c. 479 genera, and c. 6570 species. However, family and subfamily relationships in Sapindales remain unclear, making reconstruction of the order's spatio-temporal and morphological evolution difficult. In this study, we used Angiosperms353 target capture data to generate the most densely sampled phylogenetic trees of Sapindales to date, with 448 samples and c. 85% of genera represented. The percentage of paralogous loci and allele divergence was characterized across the phylogeny, which was time-calibrated using 29 rigorously assessed fossil calibrations. All families were supported as monophyletic. Two core family clades subdivide the order, the first comprising Kirkiaceae, Burseraceae, and Anacardiaceae, the second comprising Simaroubaceae, Meliaceae, and Rutaceae. Kirkiaceae is sister to Burseraceae and Anacardiaceae, and, contrary to current understanding, Simaroubaceae is sister to Meliaceae and Rutaceae. Sapindaceae is placed with Nitrariaceae and Biebersteiniaceae as sister to the core Sapindales families, but the relationships between these families remain unclear, likely due to their rapid and ancient diversification. Sapindales families emerged in rapid succession, coincident with the climatic change of the Mid-Cretaceous Hothouse event. Subfamily and tribal relationships within the major families need revision, particularly in Sapindaceae, Rutaceae and Meliaceae. Much of the difficulty in reconstructing relationships at this level may be caused by the prevalence of paralogous loci, particularly in Meliaceae and Rutaceae, that are likely indicative of ancient gene duplication events such as hybridization and polyploidization playing a role in the evolutionary history of these families. This study provides key insights into factors that may affect phylogenetic reconstructions in Sapindales across multiple scales, and provides a state-of-the-art phylogenetic framework for further research.

Inter- and intra-island speciation and their morphological and ecological correlates in Aeonium (Crassulaceae), a species-rich Macaronesian radiation.

BACKGROUND AND AIMS: The most species-rich and ecologically diverse plant radiation on the Canary Islands is the Aeonium alliance (Crassulaceae). In island radiations like this, speciation can take place either within islands or following dispersal between islands. Aiming at quantifying intra- and inter-island speciation events in the evolution of Aeonium, and exploring their consequences, we hypothesized that (1) intra-island diversification resulted in stronger ecological divergence of sister lineages, and that (2) taxa on islands with a longer history of habitation by Aeonium show stronger ecological differentiation and produce fewer natural hybrids. METHODS: We studied the biogeographical and ecological setting of diversification processes in Aeonium with a fully sampled and dated phylogeny inferred using a ddRADseq approach. Ancestral areas and biogeographical events were reconstructed in BioGeoBEARS. Eleven morphological characters and three habitat characteristics were taken into account to quantify the morphological and ecological divergence between sister lineages. A co-occurrence matrix of all Aeonium taxa is presented to assess the spatial separation of taxa on each island. KEY RESULTS: We found intra- and inter-island diversification events in almost equal numbers. In lineages that diversified within single islands, morphological and ecological divergence was more pronounced than in lineages derived from inter-island diversification, but only the difference in morphological divergence was significant. Those islands with the longest history of habitation by Aeonium had the lowest percentages of co-occurring and hybridizing taxon pairs compared with islands where Aeonium arrived later. CONCLUSIONS: Our findings illustrate the importance of both inter- and intra-island speciation, the latter of which is potentially sympatric speciation. Speciation on the same island entailed significantly higher levels of morphological divergence compared with inter-island speciation, but ecological divergence was not significantly different. Longer periods of shared island habitation resulted in the evolution of a higher degree of spatial separation and stronger reproductive barriers.

Different molecular changes underlie the same phenotypic transition: Origins and consequences of independent shifts to homostyly within species.

The repeated transition from outcrossing to selfing is a key topic in evolutionary biology. However, the molecular basis of such shifts has been rarely examined due to lack of knowledge of the genes controlling these transitions. A classic example of mating system transition is the repeated shift from heterostyly to homostyly. Occurring in 28 angiosperm families, heterostyly is characterized by the reciprocal position of male and female sexual organs in two (or three) distinct, usually self-incompatible floral morphs. Conversely, homostyly is characterized by a single, self-compatible floral morph with reduced separation of male and female organs, facilitating selfing. Here, we investigate the origins of homostyly in Primula vulgaris and its microevolutionary consequences by integrating surveys of the frequency of homostyles in natural populations, DNA sequence analyses of the gene controlling the position of female sexual organs (CYPᵀ), and microsatellite genotyping of both progeny arrays and natural populations characterized by varying frequencies of homostyles. As expected, we found that homostyles displace short-styled individuals, but long-style morphs are maintained at low frequencies within populations. We also demonstrated that homostyles repeatedly evolved from short-styled individuals in association with different types of loss-of-function mutations in CYPᵀ. Additionally, homostyly triggers a shift to selfing, promoting increased inbreeding within and genetic differentiation among populations. Our results elucidate the causes and consequences of repeated transitions to homostyly within species, and the putative mechanisms precluding its fixation in P. vulgaris. This study represents a benchmark for future analyses of losses of heterostyly in other angiosperms.

Plant growth forms dictate adaptations to the local climate.

Adaptive radiation is a significant driver of biodiversity. Primarily studied in animal systems, mechanisms that trigger adaptive radiations remain poorly understood in plants. A frequently claimed indicator of adaptive radiation in plants is growth form diversity when tied to the occupation of different habitats. However, it remains obscure whether morphological adaptations manifest as growth form diversity per se or as its constituent traits. We use the classic Aeonium radiation from the Canary Islands to ask whether adaptation across climatic space is structured by growth form evolution. Using morphological sampling with site-associated climate in a phylogenetic context, we find that growth forms dictate adaptations to the local environment. Furthermore, we demonstrate that the response of specific traits to analogous environments is antagonistic when growth forms are different. This finding suggests for the first time that growth forms represent particular ecological functions, allowing the co-occurrence of closely related species, being a product of divergent selection during evolution in sympatry.

The Implications of Incongruence between Gene Tree and Species Tree Topologies for Divergence Time Estimation.

Phylogenetic analyses are increasingly being performed with data sets that incorporate hundreds of loci. Due to incomplete lineage sorting, hybridization, and horizontal gene transfer, the gene trees for these loci may often have topologies that differ from each other and from the species tree. The effect of these topological incongruences on divergence time estimation has not been fully investigated. Using a series of simulation experiments and empirical analyses, we demonstrate that when topological incongruence between gene trees and the species tree is not accounted for, the temporal duration of branches in regions of the species tree that are affected by incongruence is underestimated, whilst the duration of other branches is considerably overestimated. This effect becomes more pronounced with higher levels of topological incongruence. We show that this pattern results from the erroneous estimation of the number of substitutions along branches in the species tree, although the effect is modulated by the assumptions inherent to divergence time estimation, such as those relating to the fossil record or among-branch-substitution-rate variation. By only analyzing loci with gene trees that are topologically congruent with the species tree, or only taking into account the branches from each gene tree that are topologically congruent with the species tree, we demonstrate that the effects of topological incongruence can be ameliorated. Nonetheless, even when topologically congruent gene trees or topologically congruent branches are selected, error in divergence time estimates remains. This stems from temporal incongruences between divergence times in species trees and divergence times in gene trees, and more importantly, the difficulty of incorporating necessary assumptions for divergence time estimation. [Divergence time estimation; gene trees; species tree; topological incongruence.].

Identifying gaps in the photographic record of the vascular plant flora of the Americas.

Field photographs of plant species are crucial for research and conservation, but the lack of a centralized database makes them difficult to locate. We surveyed 25 online databases of field photographs and found that they harboured only about 53% of the approximately 125,000 vascular plant species of the Americas. These results reflect the urgent need for a centralized database that can both integrate and complete the photographic record of the world's flora.

Prevalence and relationship of endosymbiotic Wolbachia in the butterfly genus Erebia.

BACKGROUND: Wolbachia is an endosymbiont common to most invertebrates, which can have significant evolutionary implications for its host species by acting as a barrier to gene flow. Despite the importance of Wolbachia, still little is known about its prevalence and diversification pattern among closely related host species. Wolbachia strains may phylogenetically coevolve with their hosts, unless horizontal host-switches are particularly common. We address  these issues in the genus Erebia, one of the most diverse Palearctic butterfly genera. RESULTS: We sequenced the Wolbachia genome from a strain infecting Erebia cassioides and showed that it belongs to the Wolbachia supergroup B, capable of infecting arthropods from different taxonomic orders. The prevalence of Wolbachia across 13 closely related Erebia host species based on extensive population-level genetic data revealed that multiple Wolbachia strains jointly infect all investigated taxa, but with varying prevalence. Finally, the phylogenetic relationships of Wolbachia strains are in some cases significantly associated to that of their hosts, especially among the most closely related Erebia species, demonstrating mixed evidence for phylogenetic coevolution. CONCLUSIONS: Closely related host species can be infected by closely related Wolbachia strains, evidencing some phylogenetic coevolution, but the actual pattern of infection more often reflects historical or contemporary geographic proximity among host species. Multiple processes, including survival in distinct glacial refugia, recent host shifts in sympatry, and a loss of Wolbachia during postglacial range expansion seem to have jointly shaped the complex interactions between Wolbachia evolution and the diversification of its host among our studied Erebia species.

Speciation through chromosomal fusion and fission in Lepidoptera.

Changes in chromosome numbers may strongly affect reproductive barriers, because individuals heterozygous for distinct karyotypes are typically expected to be at least partially sterile or to show reduced recombination. Therefore, several classic speciation models are based on chromosomal changes. One import mechanism generating variation in chromosome numbers is fusion and fission of existing chromosomes, which is particularly likely in species with holocentric chromosomes, i.e. chromosomes that lack a single centromere. Holocentric chromosomes evolved repeatedly across the tree of life, including in Lepidoptera. Although changes in chromosome numbers are hypothesized to be an important driver of the spectacular diversification of Lepidoptera, comparative studies across the order are lacking. We performed the first comprehensive literature survey of karyotypes for Lepidoptera species since the 1970s and tested if, and how, chromosomal variation might affect speciation. Even though a meta-analysis of karyological differences between closely related taxa did not reveal an effect on the degree of reproductive isolation, phylogenetic diversification rate analyses across the 16 best-covered genera indicated a strong, positive association of rates of chromosome number evolution and speciation. These findings suggest a macroevolutionary impact of varying chromosome numbers in Lepidoptera and likely apply to other taxonomic groups, especially to those with holocentric chromosomes. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

Phylogeny, evolution, and biogeographic history of Calandrinia (Montiaceae).

PREMISE OF THE STUDY: Calandrinia are small, succulent herbs that vary broadly in habitat, morphology, life history, and photosynthetic metabolism. The lineage is placed within the Montiaceae, which in turn is sister to the rest of the Portulacineae (Caryophyllales). Calandrinia occupy two distinct biogeographic regions, one in the Americas (~14 species), and one in Australia (~74 species). Past analyses of the Montiaceae present conflicting hypotheses for the phylogenetic placement and monophyly of Calandrinia, and to date, there has been no molecular phylogenetic analysis of the Australian species. METHODS: Using a targeted gene enrichment approach, we sequenced 297 loci from multiple gene families across the Montiaceae, including all named and 16 putative new species of Australian Calandrinia, and the enigmatic monotypic genus Rumicastrum. KEY RESULTS: All data sets and analyses reject the monophyly of Calandrinia, with Australian and New World Calandrinia each comprising distinct and well-supported clades, and Rumicastrum nested within Australian Calandrinia. We provide the first well-supported phylogeny for Australian Calandrinia, which includes all named species and several phrase-named taxa. CONCLUSIONS: This study brings much needed clarity to relationships within Montiaceae and confirms that New World and Australian Calandrinia do not form a clade. Australian Calandrinia is a longtime resident of the continent, having diverged from its sister lineage ~30 Ma, concurrent with separation of Australia from Antarctica. Most diversification occurred during the middle Miocene, with lowered speciation and/or higher extinction rates coincident with the establishment of severe aridity by the late Miocene.

Molecular evolution of key metabolic genes during transitions to C4 and CAM photosynthesis.

PREMISE OF THE STUDY: Next-generation sequencing facilitates rapid production of well-sampled phylogenies built from very large genetic data sets, which can then be subsequently exploited to examine the molecular evolution of the genes themselves. We present an evolutionary analysis of 83 gene families (19 containing carbon-concentrating mechanism (CCM) genes, 64 containing non-CCM genes) in the portullugo clade (Caryophyllales), a diverse lineage of mostly arid-adapted plants that contains multiple evolutionary origins of all known photosynthesis types in land plants (C3 , C4 , CAM, C4 -CAM, and various intermediates). METHODS: We inferred a phylogeny of 197 individuals from 167 taxa using coalescent-based approaches and individual gene family trees using maximum likelihood. Positive selection analyses were conducted on individual gene family trees with a mixed effects model of evolution (MEME). We devised new indices to compare levels of convergence and prevalence of particular residues between CCM and non-CCM genes and between species with different photosynthetic pathways. KEY RESULTS: Contrary to expectations, there were no significant differences in the levels of positive selection detected in CCM versus non-CCM genes. However, we documented a significantly higher level of convergent amino acid substitutions in CCM genes, especially in C4 taxa. CONCLUSIONS: Our analyses reveal a new suite of amino acid residues putatively important for C4 and CAM function. We discuss both the advantages and challenges of using targeted enrichment sequence data for exploratory studies of molecular evolution.

Multiple Geographical Origins of Environmental Sex Determination enhanced the diversification of Darwin's Favourite Orchids.

Environmental sex determination (ESD) - a change in sexual function during an individual life span driven by environmental cues - is an exceedingly rare sexual system among angiosperms. Because ESD can directly affect reproduction success, it could influence diversification rate as compared with lineages that have alternative reproductive systems. Here we test this hypothesis using a solid phylogenetic framework of Neotropical Catasetinae, the angiosperm lineage richest in taxa with ESD. We assess whether gains of ESD are associated with higher diversification rates compared to lineages with alternative systems while considering additional traits known to positively affect diversification rates in orchids. We found that ESD has evolved asynchronously three times during the last ~5 Myr. Lineages with ESD have consistently higher diversification rates than related lineages with other sexual systems. Habitat fragmentation due to mega-wetlands extinction, and climate instability are suggested as the driving forces for ESD evolution.

Both morph- and species-dependent asymmetries affect reproductive barriers between heterostylous species.

The interaction between floral traits and reproductive isolation is crucial to explaining the extraordinary diversity of angiosperms. Heterostyly, a complex floral polymorphism that optimizes outcrossing, evolved repeatedly and has been shown to accelerate diversification in primroses, yet its potential influence on isolating mechanisms remains unexplored. Furthermore, the relative contribution of pre- versus postmating barriers to reproductive isolation is still debated. No experimental study has yet evaluated the possible effects of heterostyly on pre- and postmating reproductive mechanisms. We quantify multiple reproductive barriers between the heterostylous Primula elatior (oxlip) and P. vulgaris (primrose), which readily hybridize when co-occurring, and test whether traits of heterostyly contribute to reproductive barriers in unique ways. We find that premating isolation is key for both species, while postmating isolation is considerable only for P. vulgaris; ecogeographic isolation is crucial for both species, while phenological, seed developmental, and hybrid sterility barriers are also important in P. vulgaris, implicating sympatrically higher gene flow into P. elatior. We document for the first time that, in addition to the aforementioned species-dependent asymmetries, morph-dependent asymmetries affect reproductive barriers between heterostylous species. Indeed, the interspecific decrease of reciprocity between high sexual organs of complementary floral morphs limits interspecific pollen transfer from anthers of short-styled flowers to stigmas of long-styled flowers, while higher reciprocity between low sexual organs favors introgression over isolation from anthers of long-styled flowers to stigmas of short-styled flowers. Finally, intramorph incompatibility persists across species boundaries, but is weakened in long-styled flowers of P. elatior, opening a possible backdoor to gene flow through intramorph pollen transfer between species. Therefore, patterns of gene flow across species boundaries are likely affected by floral morph composition of adjacent populations. To summarize, our study highlights the general importance of premating isolation and newly illustrates that both morph- and species-dependent asymmetries shape boundaries between heterostylous species.

Estimating the normal background rate of species extinction.

A key measure of humanity's global impact is by how much it has increased species extinction rates. Familiar statements are that these are 100-1000 times pre-human or background extinction levels. Estimating recent rates is straightforward, but establishing a background rate for comparison is not. Previous researchers chose an approximate benchmark of 1 extinction per million species per year (E/MSY). We explored disparate lines of evidence that suggest a substantially lower estimate. Fossil data yield direct estimates of extinction rates, but they are temporally coarse, mostly limited to marine hard-bodied taxa, and generally involve genera not species. Based on these data, typical background loss is 0.01 genera per million genera per year. Molecular phylogenies are available for more taxa and ecosystems, but it is debated whether they can be used to estimate separately speciation and extinction rates. We selected data to address known concerns and used them to determine median extinction estimates from statistical distributions of probable values for terrestrial plants and animals. We then created simulations to explore effects of violating model assumptions. Finally, we compiled estimates of diversification-the difference between speciation and extinction rates for different taxa. Median estimates of extinction rates ranged from 0.023 to 0.135 E/MSY. Simulation results suggested over- and under-estimation of extinction from individual phylogenies partially canceled each other out when large sets of phylogenies were analyzed. There was no evidence for recent and widespread pre-human overall declines in diversity. This implies that average extinction rates are less than average diversification rates. Median diversification rates were 0.05-0.2 new species per million species per year. On the basis of these results, we concluded that typical rates of background extinction may be closer to 0.1 E/MSY. Thus, current extinction rates are 1,000 times higher than natural background rates of extinction and future rates are likely to be 10,000 times higher.

A tale of two morphs: modeling pollen transfer, magic traits, and reproductive isolation in parapatry.

The evolution of the flower is commonly thought to have spurred angiosperm diversification. Similarly, particular floral traits might have promoted diversification within specific angiosperm clades. We hypothesize that traits promoting the precise positional transfer of pollen between flowers might promote diversification. In particular, precise pollen transfer might produce partial reproductive isolation that facilitates adaptive divergence between parapatric populations differing in their reproductive-organ positions. We investigate this hypothesis with an individual-based model of pollen transfer dynamics associated with heterostyly, a floral syndrome that depends on precise pollen transfer. Our model shows that precise pollen transfer can cause sexual selection leading to divergence in reproductive-organ positions between populations served by different pollinators, pleiotropically causing an increase in reproductive isolation through a "magic trait" mechanism. Furthermore, this increased reproductive isolation facilitates adaptive divergence between the populations in an unlinked, ecologically selected trait. In a different pollination scenario, however, precise pollen transfer causes a decrease in adaptive divergence by promoting asymmetric gene flow. Our results highlight the idea that magic traits are not "magic" in isolation; in particular, the effect size of magic traits in speciation depends on the external environment, and also on other traits that modify the strength of the magic trait's influence on non-random mating. Overall, we show that the evolutionary consequences of pollen transfer dynamics can depend strongly on the available pollinator fauna and on the morphological fit between flowers and pollinators. Furthermore, our results illustrate the potential importance of even weak reproductive isolating barriers in facilitating adaptive divergence.

Heterostyly accelerates diversification via reduced extinction in primroses.

The exceptional species diversity of flowering plants, exceeding that of their sister group more than 250-fold, is especially evident in floral innovations, interactions with pollinators and sexual systems. Multiple theories, emphasizing flower-pollinator interactions, genetic effects of mating systems or high evolvability, predict that floral evolution profoundly affects angiosperm diversification. However, consequences for speciation and extinction dynamics remain poorly understood. Here, we investigate trajectories of species diversification focusing on heterostyly, a remarkable floral syndrome where outcrossing is enforced via cross-compatible floral morphs differing in placement of their respective sexual organs. Heterostyly evolved at least 20 times independently in angiosperms. Using Darwin's model for heterostyly, the primrose family, we show that heterostyly accelerates species diversification via decreasing extinction rates rather than increasing speciation rates, probably owing to avoidance of the negative genetic effects of selfing. However, impact of heterostyly appears to differ over short and long evolutionary time-scales: the accelerating effect of heterostyly on lineage diversification is manifest only over long evolutionary time-scales, whereas recent losses of heterostyly may prompt ephemeral bursts of speciation. Our results suggest that temporal or clade-specific conditions may ultimately determine the net effects of specific traits on patterns of species diversification.

Small and ugly? Phylogenetic analyses of the "selfing syndrome" reveal complex evolutionary fates of monomorphic primrose flowers.

One of the most common trends in plant evolution, loss of self-incompatibility and ensuing increases in selfing, is generally assumed to be associated with a suite of phenotypic changes, notably a reduction of floral size, termed the selfing syndrome. We investigate whether floral morphological traits indeed decrease in a deterministic fashion after losses of self-incompatibility, as traditionally expected, using a phylogeny of 124 primrose species containing nine independent transitions from heterostyly (heteromorphic incompatibility) to homostyly (monomorphic self-compatibility), a classic system for evolution of selfing. We find similar overall variability of homostylous and heterostylous species, except for diminished herkogamy in homostyles. Bayesian mixed models demonstrate differences between homostylous and heterostylous species in all traits, but net effects across species are small (except herkogamy) and directionality differs among traits. Strongly drift-like evolutionary trajectories of corolla tube length and corolla diameter inferred by Ornstein-Uhlenbeck models contrast with expected deterministic trajectories toward small floral size. Lineage-specific population genetic effects associated with evolution of selfing may explain that reductions of floral size represent one of several possible outcomes of floral evolution after loss of heterostyly in primroses. Contrary to the traditional paradigm, selfing syndromes may, but do not necessarily evolve in response to increased selfing.

Evolution of multicellularity coincided with increased diversification of cyanobacteria and the Great Oxidation Event.

Cyanobacteria are among the most diverse prokaryotic phyla, with morphotypes ranging from unicellular to multicellular filamentous forms, including those able to terminally (i.e., irreversibly) differentiate in form and function. It has been suggested that cyanobacteria raised oxygen levels in the atmosphere around 2.45-2.32 billion y ago during the Great Oxidation Event (GOE), hence dramatically changing life on the planet. However, little is known about the temporal evolution of cyanobacterial lineages, and possible interplay between the origin of multicellularity, diversification of cyanobacteria, and the rise of atmospheric oxygen. We estimated divergence times of extant cyanobacterial lineages under Bayesian relaxed clocks for a dataset of 16S rRNA sequences representing the entire known diversity of this phylum. We tested whether the evolution of multicellularity overlaps with the GOE, and whether multicellularity is associated with significant shifts in diversification rates in cyanobacteria. Our results indicate an origin of cyanobacteria before the rise of atmospheric oxygen. The evolution of multicellular forms coincides with the onset of the GOE and an increase in diversification rates. These results suggest that multicellularity could have played a key role in triggering cyanobacterial evolution around the GOE.